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Abstract 



Context. The disagreement between theoretical predictions and observations for surface lithium abundance in stars is a long-standing 
problem, which indicates that the adopted physical treatment is still lacking in some points. However, thanks to the recent improve- 
ments in both models and observations, it is interesting to analyse the situation to evaluate present uncertainties. 
Aims. We present a consistent and quantitative analysis of the theoretical uncertainties affecting surface lithium abundance in the 
current generation of models. 

Methods. By means of an up-to-date and well tested evolutionary code, FRANEC, theoretical errors on surface 7 Li abundance predic- 
tions, during the pre-main sequence (pre-MS) and main sequence (MS) phases, are discussed in detail. Then, the predicted surface 7 Li 
abundance was tested against observational data for five open clusters, namely Ic 2602, a Per, Blancol, Pleiades, and Ngc 2516, and 
for four detached double-lined eclipsing binary systems. Stellar models for the aforementioned clusters were computed by adopting 
suitable chemical composition, age, and mixing length parameter for MS stars determined from the analysis of the colour-magnitude 
diagram of each cluster. We restricted our analysis to young clusters, to avoid additional uncertainty sources such as diffusion and/or 
radiative levitation efficiency. 

Results. We confirm the disagreement, within present uncertainties, between theoretical predictions and 7 Li observations for standard 
models. However, we notice that a satisfactory agreement with observations for 7 Li abundance in both young open clusters and binary 
systems can be achieved if a lower convection efficiency is adopted during the pre-MS phase with respect to the MS one. 

Key words. Stars: abundances - evolution - interiors - low-mass - pre-main sequence - Hertzsprung-Russell and C-M diagrams 



1. Introduction 

In the last two decades, a large number of 7 Li observations 
have been collected for isolated stars, binary systems, and open 
clusters from the pre-M S to the late MS phases (see e.g. Table 1 
and references therein in Jeffries 2000; Sestito & Randich 2005), 
showing that 7 Li depletion is a strong function of both mass 
and ag e. A detailed and homogen eous analysis has been carried 
out by ISestito"& Randich (2005), who determined surface 7 Li 
abundance for a large sample of open clusters in a wide range of 
ages and chemical compositions, supplying a useful tool for ac- 
curately analysing the temporal evolution of surface 7 Li abund- 
ance. 

Open clusters and detached double-lined eclipsing binaries 
(EBs) are ideal systems for testing the validity of stellar evol- 
utionary models, since their members have the same chemical 
composition and age. As a consequence, they allow the differ- 
ent lithium depletion pattern to be investigated as a function of 
the stellar mass once the age and the chemical composition have 
been kept fixed. 

Besides the large amount of 7 Li data available, a strong 
effort in theoretical modelling has been made in the past years, 
and many different theoretical scenarios have been proposed 
to explain the observed surface 7 Li abu ndance and its tem- 
poral evolution (see e.g. the reviews in Delivan nis et al.l f2000: 



iPinsonneault et alJ l2000t ICharbonnel et alJ l2000l) . both in the 



Pinsonneault et al. 


1990; Pinsonneault 


1994; Chabover et al. 


1995; D'Antona & Mazzitelli 1997; 


Ventura et al. 


1998; 


Piau & Turck-Chieze 2002; D'Antona & Montalban 


2003; 


Montalban & D'Antona 2006). 



Standard models assume a spherically symmetric structure 
and convection and diffusion are the only processes that mix sur- 
face elements with the interior. Although the validity of such 
models in reproducing the main evolutionary parameters has 
been largely tested against observations, they fail to reproduce 
the observed 7 Li abundances. Indeed, standard models show a 
7 Li depletion during the pre-MS phase that is much stronger 
than observed, wh ile the opposite occurs in the MS phase (see 
e.g., [Jeffries 2000). Moreover, they cannot fully account for the 
formation of the so-called lithium dip fo r MS stars in the tem- 
perat ure range 6000 K £ ^ 7000 K dBoesgaard & Tripiccol 
119861) . see e.g. lRicher & Michaudldl993l) . 

The comparison between theory and observation is im- 
proved, in some cases, by introducing non-standard processes 
into the models, e.g. rotation, gravity waves, magnetic fields, and 



accretion/mass loss (Pinsonneault et al. 

I ' — 1 1 rr^^^^^^^^^™ 

Chabover et al. 1995; Talon & Charbonnel 
1998 1 : iMendes et all 119991: ISiess et all Il999t 
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1990; D'Antona 1993: 
" 19981: 1 Ventura et al 
D'Antona et al 



2000; CharbonML&sIo3 12005; Baraffe & Chabrier 2010; 
Vick et all l2010h . All these processes produce structural 
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changes, with a related strong effect on lithium abundance 
(see e.g. the review s by Charb onnel etaD 120001: iTalonl 120081: 
Talon & Charbonnel 20 lTJ. In particular, models with rotation- 
induced mixing plus gravity waves are able to reproduce 7 Li 
the depletion during the MS and post MS phases (i.e. the 
lithium dip feature and red-giant branch a bundances, see e.g., 
iTalon & CharbonnefeO 10b [Pace et al.ll20l2h . 

A crucial point in stellar modelling, both for standard 
and non-standard models, concerns the treatment of the over- 
adiabatic convection efficiency in the stellar envelope, which is 
an important issue for lithium depletion, too. In evolutionary 
codes, the most widely used convection treatm ent is the sim- 
plified scheme of the mixing length theory (MLT. iBohm- Vitensel 
119581) . In this formalism, convection efficiency depends on a free 
parameter to be calibrated. It is a common approach to calib- 
rate it by reproducing the solar radius. This choice usually gives 
good agreement between models and photometric data; how- 
ever, to reproduce the effective temperature of stars with dif- 
ferent masses in different evolutionary phases, an ad hoc value 
of the mixing length para meter should be adopted, as suggested 
by observations (see e.g..|Chieffi et al.lll995b iMorel et aLlfe OOO ; 
Ferra ro et al.ll2006HYildizll2007l:lGennaro et al.ll201 UlPiau et all 
201 1; Bonaca et al. 2012) a nd detailed hydrodynam ical simula- 



tions (see e.g.,|Ludwig et ail 19991: [Trampedachl2007l) . 

The main goal of this paper is to re-examine the old lith- 
ium problem in light of the improvements in the adopted phys- 
ical inputs and observational data and to perform a quantitat- 
ive analysis of the uncertainties affecting surface lithium de- 
pletion during the pre-MS phase. The aim is to compute, by 
means of updated models, theoretical error bars to be applied 
to the comparison between predictions and data available for 



stars in 



young open cluster and bi nary systems, as partially done 

ai- *-vfrVi £»■»- \T7/-Vf-L*-c ( oao <n rr H' Ant/inn JPr A /Tot 71 1 il I 1 Q Q AV 



in earlier other works (see e.g 
Swenson et al.lll994t IVentura et al. 
2002tlSestitoetal.ll2006h 



D'Antona & Mazzitelli 1984 
Piau & Turck-Chieze 



1998; 



The paper is structured in the following way. Section [2] 
presents the adopted 7 Li data sample for the selected open 
clusters, followed by a brief description of present models (Sect. 
|3). In Sect. |4]we evaluate the main theoretical uncertainties af- 
fecting surface lithium abundance. Finally, in Sect. [5] the com- 
parison between predicted and observed lithium abundances for 
both young open clusters and binary systems is discussed. 



2. Lithium data 

Surface 7 Li abundances for young open clus ters are taken from 
the ho mogeneous database made available bv lSestito & Randichl 
(12001 . Here, we focus our analysis on clusters younger than 
about 150 - 200 Myr, in order to avoid MS depletion effects 
(see e.g.. ISestito & Randichl 120051) . with different metallicities 
for which a significant number of data in a wide range of ef- 
fective temperatures are available. The clusters that satisfy these 
criteria are, Ic 2602, a Per, Blanco 1, Pleiades, and Ngc 2516. 

Lithium abundances for young d ouble-lined eclipsing binar - 
ies are not present in the database bv lSestito & Randichl (12005 ). 
but they have been measured by different authors, as we discuss 
in Sect. 



3. Theoretical stellar models 

Present stellar models were compu ted with an updated versio n 
of the FRANEC evolutionary code dDeglTnnocenti et al.l 120081) . 
which adopts the most recent input physics, as described in de- 



tail by Tognell i et al ] (120111) . The initial deuterium mass frac- 
tion abundance is fixed to Xd = 2 X 10 as a representative 



value for population I stars (see e.g. Geiss & Gloecklerl [l 998 ; 
iLinskv et al.l 120061 Steigman et al.ll2007 ). The logarithmic ini- 
tial lithium abundance is assumed to be = 3.2 + 0.2 (see e.g., 
Jeffries 20061 lLodderseTaT] l2009). which approximatively cor- 
responds to Xiu * 7 x 10~ 9 - 1 x 10~ 8 in dependence on the 
metallicity adopted for the modelfl 

Convection is treated according to the mixing length th eory, 
using the same formalism presented in lCox & Giulil 11968). The 
adopted reference value of mixing length parameter is a — 1.0 
(as suggested by present comparison with pre-MS data, see Sect. 
ED. 

4. Theoretical uncertainties 

4. 1 . Chemical composition 

To properly calculate pre-MS evolution suitable initial abund- 
ances of helium, light elements, and metals are needed. For most 
of the stars, however, only the [Fe/H] value is available, so the- 
oretical or semi-empirical assumptions are required. 

Assuming for Populatio n I stars a solar-sca led heavy ele- 
ments distribution (see e.g., Asp lund et al.ll2009l for a detailed 
review), the Z/X value currently present at the stellar surface 
can be directly inferred from the observed [Fe/H]. For all the 
stars analysed in the paper, this value can be safely adopted as 
a good approximation of the initial one over the whole struc- 
ture, since the effect of microscopic diffusion is negligible ow- 
ing to the very young ages involved. The initial helium content 
of the star Y cannot be directly measured in the stellar spectra 
of cool stars, so a further relation between the initial metallicity 
and helium of the star is required. A common wa y to proceed is 
to ass ume the following linear relation (see e.g., Gennaro et al. 
I2010h : 



AY 

Y = Y P + — Z 
AZ 



(1) 



where Yp and AF/AZ represent, respectively, the primordial he- 
lium abundance and the helium-to-metal enrichm ent ratio. For 
the ca lculations we adopt Y „ = 0.2485 + 0.0008 (ICvburt et al.l 
12004 and AF/AZ = 2 + 1 Casagrande et al.ll2007l) . Thus, the 
metallicity of the star can be obtained directly from the follow- 
ing equation, 



(1 - F P )(Z/X) e 



10-[Fe/H] _ (1 + AY/AZ)(Z/X) e 



(2) 



once the solar (Z/X) G has been specified. Regarding this 
last quantity, there are several values a dopted b y differ- 
ent au thors, i.e. the still widely adopted Grevesse & Noels 
(fl99l (GN93, (Z/X) B = 0.0244), iGrevesse & Sauvall (119981) 
(GS 98,(Z/X)o = 0.0 231) and the recent determinations 
bv lAsplundetalJ d2005b (AS05, (Z/X) Q = 0.0165) and 
Aspl und et al.l (l2009f) (AS09. (Z/X) e = 0.0181), which are based 
on detailed 3D hy drodynamical atmosphere models. Recently, 
ICaffauetal.l d2010) (CLIO) have found a value for the solar car- 
bon photospheric abundance higher by about 0. 1 dex than the 
previous one derived by AS09. This also leads to an increase 
in the solar metallicity-to-hydrogen ratio, namely (Z/X) Q = 



1 We adopt a simple scaling of the initial 7 Li abundance with the 
metallicity because we are mainly interested in reproducing the lithium 
depletion pattern, i.e. e Li - T eff , which is independent of the initial 7 Li 
abundance. 
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Figure 1. Relative differences of the Rosseland radiative opacity coefficients (A/cr/kr) computed by adopting the AS05 and GN93 
(left panel), AS05 and GS98 (central panel), and AS05 and AS09 (right panel) solar mixtures, for Z =0.0129 and Y = 0.274. The 
three panels also show the region covered by the entire convective envelope of masses in the range 0.6 - 1.2 M (thick-dashed box), 
for the same chemical composition of the opacity tables and adopting the mixing length value a = 1 .00. 



0.0211, which is higher than the AS09 and much closer to the 
GS98 one. 

Our models are computed adopting the AS05 mixture, for 
consistency with the extended pre-MS tra cks and isochrone s 
database already provided by our grou p0 (Togn elli etal]l20TTh . 
whereas for the conversion of [Fe/H] into (Y, Z) we prefer to 
use (Z/X) Q = 0.0181 from the more recent AS09 heavy ele- 
ments distribution. The inconsistency that may arise is negli- 
gible. Indeed, we verified that the effect on pre-MS models of 
adopting the AS05 or AS09 distribution in the opacity, once Z 
and F have been kept fixed, is much lower than the variation pro- 
duced by a change of (Y, Z) related to the error on (Z/X) (see 
iTognelli et al.ll201 H and the following discussion). 

For the uncertainty on (Z/X) e , the commonly sugges- 
ted value is about ±15 % (see e.g., iBahcall et ail l2004t 
Bahcall & Sereneliil l2005h . However, to take the difference 
between recent (Z/X) Q determinations into account, i.e. GS98, 
AS09, and CLIO, a larger uncertainty of about +25% with re- 
spect to AS09 is required. Thus, we use a final uncertainty of 
+25/-15%on(Z/X) . 

Besides the uncertainties on Yp, AF/AZ, and (Z/X) Q , the ini- 
tial F and Z abundances are obviously also affected by the ob- 
servational error on [Fe/H]. Generally the errors quoted in the 
literature vary from about +0.01, probably underestimated, to 
±0.1. Here we adopt as a conservative error A[Fe/H] = ±0.05. 

By means of eqs. ([TJ and ©, taking the quoted uncertainties 
into account, we obtain eight values of (F, Z) for each of the 
selected clusters, for which we compute pre-MS models. More 
precisely, the different F and Z values are calculated by adopting, 
in turn, the minimum and the maximum values of one of the four 
parameters ([Fe/H], F P , AF/AZ, and (Z/X) Q ), while the others 
are kept fixed to their central value. We computed two additional 
models with the maximum and minimum values of the estimated 
initial 7 Li abundance, which, as already mentioned in Sect. [3] is 
set to eu - 3.2 ± 0.2. Obviously a change in the initial chemical 

2 The database contains a very large grid of pre-MS mod- 
els and isochrones between 1-100 Myr (for more details see, 
Tognelli et al. 2011). The corresponding database is available at: 
http : //astro . df . unipi . it/stellar-models/ 



composition also affects the position of the star in the colour- 
magnitude diagram, hence the age and mixing length parameter 
determination. This effect has been taken into account (see Sect. 

ED- 

Another source of uncertainty related to the stellar chemical 
composition that has to be considered is the assumed distribu- 
tion of hea vy elements, at fix ed Z, which strongly affects opa- 
city values (ISestito et al.ll2006T) . Since the opacity determines the 
temperature gradients and thus the extension of the convective 
envelope (in mass and temperature), a variation in this quant- 
ity, due to the uncertainty on the adopted mixture, modifies the 
lithium-burning rate and its resulting surface abundance. 

Figure Q] shows the relative differences among the OPAL ra- 
diative opacity tabled computed by adopting the GN93, GS98, 
AS05 (our reference model), and AS09 solar mixtures. To make 
the figure much clearer , we also show a box representing the re- 
gion covered by the entire convective envelope for stellar models 
with masses in the range 0.6 - 1.2 M and [Fe/H] = +0.0. The 
model ages have been chosen in the range 1-10 Myr, which 
roughly corresponds to the phase of efficient lithium burning for 
such masses. 

The upper panel of Fig. [2] shows the change in surface lith- 
ium abundance resulting from the adoption of the aforemen- 
tioned solar mixtures in the opacity tables, in stars with differ- 
ent masses, namely 0.6, 0.8, 1.0, and 1.2 M . The increase of 
the iron group elements at fixed metallicity Z leads to higher 
radiative opacity and in turn to a deeper convective envelope 
and consequently to greater lithium depletion. This is exactly 
what occurs when updating the heavy elements mixture to the 
recent AS05 or AS09 from the older GN93 or GS98. The higher 
the mass, the lower is the surface lithium depletion and, con- 
sequently, the sensitivity of the lithium burning to the opacity 
change. Since AS05 and AS09 metals distributions are quite 



3 We use the OPAL opacity table released in 2005 
for logT[K] > 4.5, which are available at the url, 
http://opalopacity.llnl.gov/opal.html For lower tem- 
peratures we use the Ferguson et al. (2005) radiative opacity, available 
at the url, http : //webs .wichita. edu/physics/opacity/ 
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Figure 2. Comparison among the surface lithium abundance ob- 
tained with our reference set of tracks (solid line) and models 
with different assumptions on the adopted physical inputs, for 
M = 0.6, 0.8, 1.0, and 1.2 M with Z = 0.01291, Y = 0.274, 
and a = 1. Upper panel: effect of the change of the solar mix- 
ture (GN93, GS98, AS05, and AS09 ones) in the opacity tables. 
Bottom panel: effect of adopting the OPAL 2006 and PTEH EOS. 



similar, we expect negligible differences in the opacity coeffi- 
cients, and thus in the model predictions. 

4.2. Opacity coefficients. 

Besides the influence on the opacity of the heavy element dis- 
tribution, it is worth analysing the error on the calculation 
of Rosseland radiative opacity coefficients kr at fixed chem- 
ical composition. The current version of the opacity tables 
we adopt in present calculations (i.e. OPAL 2 005 see e.g., 
Rogers & Igles ial[l992l : llglesias & Ro gers 1996) does not con- 
tain information about the related uncertainty. Thus to give a 
conservative uncertainty estimation on kr, we evaluated the rel- 
ative diff erences between the OPAL and the OP ( Opacity Project 
see e.g. ISeaton et alJ 1199-41: iBadnell et al J 12005) radiative opa- 
city coefficients in their full range of validity, once the same 
chemical composition has been adopted. We found that the 
maximum/minimum relative difference between the two opa- 
city tables is close to ±5% in the region of interest for the 
present calculations (i.e. the convective envelope) (see also, 
Neuforge-Verheecke et"aIll200lHBadnell et al-lkooa I Valle et alJ 
20121) . Thus, we assume the value of Akr/a-r = +5% as a con- 
servative uncertainty. 



4.3. Equation of state 



Owing to the complexity of the evaluations of the various ther- 
modynamical quantities, which are strictly correlated among 
each other, it is very difficult to assess a precise uncertainty on 
the EOS tables. An idea of how the current indetermination on 
the EOS propagates into stellar evolutionary predictions can be 
obtained by computing models with two different EOS tables 
that have been widely adopte d, namely the O PAL EOS 2006 
(our reference one) and PTEH dPols et al.| [T995fl The compar- 
ison between the OPAL and the PTEH is useful to assess the ef- 
fect of adopting a completely different treatment of the gas in 
stellar conditions, the two EOS being computed, respectively, 
in the formalism of the physical and chemical picture (see e.g., 
iTrampedach et al.ll2006l) . 

The influence of the adopted EOS on the location in the 
HR diagram has already been dis cussed in several papers for 
both pre-MS evolut i on (se e e.g., iMazzitellil 119891: iD'Antonal 
1993; Tognell i et al.l 120 111) and low-mass MS stars (see e.g. , 
Dor manetalJ 119891: iNeecd 11984 IChabrier & Baraffd 119971: 
di Criscienzo et alj|2010t) . Here, we simply recall that the mod- 
els are particularly affected by the EOS in all the phases where 
a thick convective envelope is present, i.e. the pre-MS or MS 
structures of low- and very low-mass stars. In these phases, when 
lithium burning is efficient, the resulting surface 7 Li abundance 
is quite sensitive to the adopted EOS, too, as shown in the bottom 
panel of Fig. [2] As noticed above, surface lithium abundance gets 
less affected by the EOS change as the mass increases, because 
of the progressively reduction of the burning efficiency. 



4.4. 7 Li(p,a)a cross section 



Lithium destruction is obviously dependent on the 7 Li(p,a)a 
cross section. However, the current uncertainty on the quoted 
react ion rate for bare nuclei is quite small ( a few percent see 
e.g.. lAngulo et alJ[l999T: lLattuada et al.ll2001h . so that the effect 
on 7 Li abundance of such error is very small compared to the 
other error sources. We adopt the value of +5% as a conservat- 
ive uncertainty on this quantity. 



4.5. Total uncertainty on 7 Li surface abundance predictions 



The partial uncertainty due to each parameter/physical input was 
obtained by the difference between the reference model, which 
is the one computed with the reference values of all the paramet- 
ers, and the model computed by varying such parameter. This 
procedure was iterated for all the uncertainty sources discussed 
in the text. Then, the total error on surface lithium abundance 
predictions was computed by quadratically adding all the par- 
tial errors. We want to emphasize that the uncertainty analysis 
was performed for all the chemical compositions suitable for 
the selected clusters. Thus, for each cluster, error bars consist- 
ent with its chemical composition, mixing length parameter, and 
age were evaluated. 



4 



E. Tognelli et al.: 7 Li surface abundance in pre-MS stars. 



Table 1. Main properties adopted/derived for the five selected 
open clusters. 



Cluster: 



[Fe/H]: 



(Y, Z): 



age (Myr): age (Myr): 
(l ov = 0.0) (A ov = 0.2) 



Ic2602 


+0.00" 


(0.274, 0.0129) 


1.68 ±0.1 


40 ± 10 


55 ± 10 


a Per 


-0.10'' 


(0.269, 0.0104) 


1.68 ±0.1 


60 ± 10 


75 ± 10 


Blanco 1 


+0.04 c 


(0.276, 0.0141) 


ISO** 


110 ±30 


130 ±30 


Pleiades 


+0.03'' 


(0.276, 0.0138) 


1.90^;i 


120 ±20 


130 ±20 


Ngc2516 


-0.10 s 
+0.07 f 


(0.269, 0.0104) 
(0.278, 0.0150) 


1.90 ±0.1 
1.90 ±0.1 


130 ±20 
130 ± 20 


145 ± 20 
145 ± 20 



Notes. The columns list, respectively, the cluster's name, [Fe/H], initial 
helium and metal abundance (Y, Z), the mixing length parameter calib- 
rated on MS stars (a-Ms), the best fit age without core overshooting, and 
the best fit age with a core overshooting parameter set to A ov = 0.2. 



Reference s. (fl) iD'Orazi & RandicU i2009h: { b) Balachandra n et al 
20TTh : (c) lFord et all <2005h: w ISoderblom etal] ( 120091) : M Sung ct al 



2002); (/) Terndr up et alJ J2002h . 



5. Surface lithium abundance: theory vs 
observations 

5. 1 . Young open clusters 

The clusters age and the mixing length parameters for MS 
stars are determined by comparing the observed CMDs with 
the present theoretical isochrones. The age is largely affected 
by the lack of stars near the overall contraction region of such 
young clusters, and it is only marginally affected by the uncer- 
tainty on the chemical composition adopted for the calculations. 
Similarly, the uncertainty on the calibration of qtms essentially 
comes from the spread of the MS in the CMD. The effects of the 
indetermination on age and cms on surface 7 Li abundance are 
evaluated for each cluster and quadratic ally added to the other 
error sources to define the theoretical surface lithium abundance 
error bars. Table [TJ summarizes the results for the main para- 
meters of each cluster, namely, chemical composition, otms, and 
age. For Ngc 2516 two different [Fe/H] values are available: 
(1) [Fe/H] = -0.10 (photometric) and (2) [Fe/H] = +0.07 (spec- 
troscopic). Since the spectroscopic value is still quite uncertain, 
because it has been determine d using only two stars (see the dis- 
cussion in Terndru p et al.l200"2l) . we decided to also use the pho- 
tometric one for the following comparison with surface lithium 
data. 

The age determination for the selected clusters was per- 
formed both with models without core overshooting (A m = 0) 
and with a standard upp er limit for the presence of overs hoot- 
ing, A m = 0.2 (see e.g.. iBrocato et al]l2003t IClaretl I2007L and 
references therein), in order to take a suitable range of A m into 
account. However, as shown in Table[TJ the age difference is for 
all the cases close to the quoted uncertainty. Moreover, the im- 
pact of the age uncertainty on €\±—T e s profile is almost negligible 
for clusters older than about 80 - 90 Myr, since all the stars have 
already reached the ZAMS, so that their position in the CMD is 
weakly dependent on the age. The situation is different for very 
young clusters, such as Ic2602 and a Per. In this case low-mass 



4 The PTEH EOS was computed using the FreeEOS fortran library 
developed by A.W. Irwin, which allows computing the EOS by the free- 
energy minimization technique. One of its advantages is the possibility 
of setting several flags to mimic other historical EOS. 



stars are close but not yet in ZAMS, so an age variation produces 
an appreciable change in T e g. 

As for the age, the ctms uncertainty also has a small effect on 
the predictions of surface 7 Li, since in the age and mass range we 
are investigating, the stars undergo an almost negligible deple- 
tion during MS. Only the surface 7 Li abundances for the lowest 
masses of the sample, i.e. M SS 0.6 M Q , are marginally affected by 
such a variation in the predicted 7 Li. Besides this, q;ms slightly 
changes the effective temperature of the stars, which, however, is 
a second-order effect when compared to the r er r shift introduced 
by the other uncertainties previously discussed. 

Regarding the comparison between surface 7 Li data and 
present models, as a first step, we make the assumption of a con- 
stant value of the mixing length parameter from the early pre- 
MS phase to the MS. Figure [3] shows the comparison between 
the theoretical predictions and data for each cluster. As one can 
see, the models with ctpms = <*ms (dashed line with filled black 
squares) fail in reproducing the observed 7 Li abundances in al- 
most all the selected clusters, for stars less massive than about 
1 Mq, even if theoretical and observational errors are taken into 
account. For these stars, the predicted 7 Li is systematically lower 
than what is observed, with differences as great as 1 dex for 
low-mass stars (about 0.6 - 0.7 M ), confirming the well known 
disagreement between theory and observations for 7 Li surface 
abundance, even in light of the recent data and taking the updated 
theoretical models into account within present error estimates. 

Models computed with czpms = <*ms partially agree with data 
only in the case of a Per for M <: 0.7 M G , and Ngc 2516 if the 
low photometric [Fe/H] value is adopted (bottom left panel of 
Fig. [3}. If the spectroscopic [Fe/H] value is used for NGC 2516, 
the predictions, as the other clusters, do not match the observa- 
tions for M £ 1 M . However, we emphasize that, for these two 
clusters, models and data are compatible each other because of 
the large 7 Li abundance scatter present among stars with similar 
T e ff (about 1 dex), combined with the large error bars on theor- 
etical predictions. 

Since in most of the cases the models with q;pms = c^ms dis- 
agree with the data, and given the high sensitivity of 7 Li surface 
abundance predictions to the convection efficiency, it is worth 
exploring the possibility that the mixing length parameter value 
varies from the pre-MS to the MS phases. Indeed, a possible de- 
pendence of a on the evolutionary phase (and/or gravity, T e s, 
mass) is suggested from both observations and hydrodynamical 
simulations, as discussed in the introduction. Thus, we computed 
models with different values of apMS, namely, ccpms = 1-0, 1.2, 
1.4, 1.68, and 1.9, once q-ms and the ages have been fixed by 
the comparison in the CMD. Figure [3] shows the comparison 
between our 'best fit' models and 7 Li data for each cluster (dot- 
ted lines and filled red squares). The theoretical error bars com- 
puted for each cluster are also shown. 

We emphasize that a satisfactory agreement with all the 
clusters in the sample (with the exception of the Pleiades) can be 
achieved by assuming the same pre-MS convection efficiency, 
namely apMS = 1 -0. Such low-convection efficiency models are 
able to reproduce, within the error bars, the mean depletion pro- 
file even for low-mass stars, especially in the case of Ic 2602 and 
a Per. 

As shown in Fig. [3] the poorest match between theory 
and data is achieved for the Pleiades. The hottest stars are 
nearly compatible within the error bars with the observations, 
which show a surface abundance about 0.2 - 0.3 dex lower 
than the predicted one. A possible way to improve the agree- 
ment with these stars is to adopt an initial lithium abund- 
ance of about fLi ~ 3. However, this method does not im- 
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Figure 3. Comparison between our model predictions and observational data (filled circles) for surface lithium abundance, for the 
selected sample of young open clusters, namely (from the top left panel), Ic2602, a Per, Blanco 1, Pleiades, and Ngc 2516. In 
each panel we show both the low-convection efficiency (o^ms = 1 -0, red-filled squares and dotted line) and the high-convection 
efficiency models (ffpMS = i^ms, dashed line and small black-filled squares). We also plotted the mass and the corresponding 
theoretical uncertainties on both the effective temperature and lithium abundance, for low- and high-convection efficiency models. 



prove the agreement with the low-mass star s, a problem still 
largely discus s ed in the litera t ure (s e e e.g. King et al. feOQO; 



Jeffri es! l200d lUmezu & Said 120001: iD'Anton a & Montalban 
2003HClarke et aUl2004HXiong & Dengll2006HKing et alJl2010l 



and references therein). 



The results we obtain for young open clusters confirm the 
partial results of previous analysis, which have noticed that 
models with low-convection efficiency during pre-MS phase 
agree much better with lithium observati on than those with 
solar or MS calibrated values (see e.g.. I Ventura et al] [l998; 
ID'Antona & Mont alban 2003t lLandin et al.ll2006l) . 
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Figure 4. HR diagram for the selected sample of EBs (see text). For each system we plotted the models for the labelled chemical 
composition and mass, corresponding to the primary (blue line) and the secondary star (red line). Models have been computed 
by adopting three different values for the mixing length parameter, namely q-pms = 1.0 (solid line), 1.20 (dashed line), and 1.68 
(dot-dashed line). Observational data are shown (filled circles), together with their uncertainty. 



5.2. Binary stars 

Binary systems and, in particular, the subclass of detached 
double-lined eclipsing binaries (EBs) are severe tests for stel- 
lar models. Indeed, for EBs independent measurements of mass, 
radius, and ef fective temperature a re available (for a detailed re- 
view see e.g jMathieu et al.1l2007l) . 

The validity of our theoretical models have already been 
tes ted against a l a rge sa mple of pre-MS binaries (26 objects) 
by iGennaro et al.l d201 ll) . using the models of the Pisa pre-MS 
database against observations by means of a Bayesian method. 
The present pre-MS models differ from those available in the 
quoted database only in the minimum value of the mixing length 
parameter, i.e. a — 1.0 instead of a — 1.2. 

Fro m the sample of EBs presented in IGennaro et al.l 
d20Tlh . we selected a subsample of binary systems for 
which surface lithium abundances are available, namely 
ASAS J05282 1+0338.5 (IStempels et al.ll2008h. EK Cep dPoppeil 
fl987b. RXJ 0529.4+ 0041 A dCovino et alj2004l) . and VI 174 Ori 
dStassun et al.l l2004). Table [2] summarizes the main parameters 
of each system: mass, effective temperature, luminosity, lithium 
abundance, and [Fe/H]. The corresponding models are computed 
forapMS = 1-0, 1.2, and 1.68. 

Figure |4] shows the HR diagram of th e four selected systems 
compared with our evolutionary tracks. IGennaro et al ] d2011h 



have already shown that theoretical models with low initial he- 
lium abundance and mixing length parameter agree better with 
the data of pre-MS binary systems, in particular for those ones 
with at least one component near the Hayashi track. In our 
sample only EK Cep does not have stars near the Hayashi track. 

For ASAS J052821+0338.5 (Fig. g), both the lowest and 
the highest otpms values are compatible with the primary star, 
whereas apMS = 1.0 - 1.2 is required to match the secondary. 
For Ek Cep we cannot constrain the mixing length value dur- 
ing the pre-MS since both stars are approaching the ZAMS, and 
consequently their position in the HR diagram is not sensitive 
to the choice of q;pms- Moreover, we can not achieve a sat- 
isfactory agreement betwe en our model and the primary star, 
as already pointed out by IGennaro et al.l (1201 ll) . Similarly to 
ASAS J052821+0338.5, RXJ 0529.4+0041 A has two stars near 
the 'heel'. As shown in Fig. [4] the three different apMS are all 
compatible with both stellar components within the observa- 
tional uncertainties, which are quite large. 

VI 174 Or i is m uch more problematic. As discussed in 
Gennaro et al 1 (1201 lh . the two stars show a peculiar position in 
the HR diagram. None of the present models (or other models 
widely adopted in the literature) can reproduce the correct posi- 
tion of the secondary by adopting the measured mass and chem- 
ical composition. A possible explanation of such a peculiar po- 
sition in the HR diagram can be the presence of a large system- 



7 



E. Tognelli et al.: 7 Li surface abundance in pre-MS stars. 
System: ASAS J052821 .0338.5 System: EK Cep 



3.6 




1 1 1 1 1 1 1 


' I 




3.4 










3.2 




- 




3.0 








2.8 




Primary: M= 1 .387 M Q 






2.6 


hi -i 


Secondary: M=1.331 M Q 
Primary track, a PMS =1.0 






2.4 




Secondary track, a PMS =1.0 
Primary track, a PUS =« us 






2.2 




Secondory trock, a PMS =a MS 








□ 


Primary track uncertainty region, 






2.0 
1.8 


- \^\ Secondory trock uncertainty regio 

[Fe/H] = -0.15, a MS =1.68 
I I 


i a P Ms =1 -° 

. I 





5000 
T.„ [K] 









System: 


RXJ 0529.4 + 0041 A 




3.6 




i i 






3.4 


















3.2 






1 J i J 




3.0 








+ 














2.8 


: 


Primary: M= 1 .270 M 






\ 

z" 


2.6 


h I-i 


Secondary: M = 0.930 M s 






cn 

o 






Primary track, a PUS — 1-0 


t 




2.4 




Secondary track, a PMS =1.0 












Primary track, a PUS = a MS 








2.2 




Secondory trock, a PUS =a MS 










□ 


Primary track uncertainty regio 


1, a PUS =1.0 






2.0 


- [ZZ] Secondory trock uncertointy region, a PWS =1.0 






1.8 


[Fe/H] = -0.01, a MS =1.68 
I I 








5000 





3.6 




3.4 




3 2 




3.0 








2.8 


\ 




~£ 


2.6 






O 






2.4 




2.2 




2.0 



f-f-i Secondary: M=1.124 W Q 

-- Secondary track, a PUS — 1-0 
-- Secondary track, a PMS =a MS 

I | Secondory trock uncertointy region, ap M s = 1.0 
[Fe/H] = + 0.07, « us =1.68 



) 5000 450 

T eff [K] 
System: VI 1 74 Or! 





3.6 




3.4 




3.2 




3.0 


+ 






2.8 


\ 




~£ 


2.6 






o 






2.4 




2.2 




2.0 



i i i i i 


















wjH Primary: M= 1.009 M a \ 


i 
i 


: 


f-jjH Secondary: M = 0.731 W Q \ 






— Primary track, a PUS = 1.0 






-- Secondary track, n PMS =1.0 






— Primary track, a PUS =a us 






-- Secondory trock, a PUS =a MS 






\^\ Primary track uncertainty region, Q PUS =1-0 






- \^\ Secondory trock uncertointy region, a PMS =1.0 






[Fe/H] = -0.01, a us =1.68 
I I I 


. ■ ■ ■ I I.. 





5000 
T,„ [K] 



Figure 5. Comparison between theoretical predictions and observations for surface lithium abundance in the selected binary systems. 
Lithium evolutionary tracks for low-convection efficiency (ctpms = 1.0) corresponding to the primary (blue-solid line) and the 
secondary star (red-dashed line) have been computed for the labelled chemical composition and mass. Shaded area represents the 
uncertainty regions for the primary (dark-shaded area) and secondary stars in the system (light-shaded area). We also show the 
primary and secondary evolutionary tracks computed with <?pms = Q'ms (green-solid and black-dashed lines, respectively). 



atic uncertainty introduced b y the adopted spectral type-effective 
temperature scale (see e.g., iLuhrnan et al. 1997t IStassun et all 
2004; iHillenbrand & White! l2004t iGennaro et al.l |201j. To be 
close to our coolest model (q-pms = 1 -0), an increase of about 
300 K in the secondary effective temperature would be required, 
which would correspond to a primary effective temperature in- 
crement of about 400 K. However, it seems unlikely that such a 
large shift could be caused uniquely by the adoption of a inad- 
equate spectral type-effective temperature scale. 

Figure|5]shows the comparisons between theoretical and ob- 
served lithium surface abundances. The evolutionary track of 
surface lithium abundance is shown. In the case of EK Cep we 
do not show the primary because the lithium abundance is not 
currently available. The figure shows the tracks computed with 
the low-convection efficiency, i.e. q-pms = 1.0 and the models 
with apMS = a M s = 1-68. 

The lithium predictions of our models computed with both 
Q'pms = 1.0 and c^pms = ffMS = 1-68 are in good agreement 
with data for the primary components (1.0 M ^ Mj ^ 1.4 M ), 
within the uncertainties, since lithium depletion is almost neg- 
ligible for such masses. Therefore, the primary components be- 
longing to our sample do not allow further constraints on the 
Q'pms value. In contrast, the impact of otpms gets stronger and 
stronger as the mass decreases below about 1 M . The second- 



ary of EK Cep and RXJ 0529.4+0041 A might thus give use- 
ful constraints^. Unfortunately, among the selected systems, the 
lithium data for the secondary components are quite uncertain, 
with errors as large as 0.5 dex (see Table avoiding a robust 
comparison. The model with a low-convection efficiency is fully 
compatible with the data for the secondary component of EK 
Cep, although we cannot exclude ctpms = 1.68, while in the 
case of RXJ 0529.4+0041 A, nothing can be concluded owing 
to extremely uncertain lithium abundance determinat ion, which 
is on l y a lower hmit (see also the discussion in. [Alcala et al 



only 

200rt (Covinoe7aIll2001t Si Antona & Montalban 2003). From 



this analysis it is evident that an effort to improve the measure- 
ments of the main parameters of binary systems, lithium abund- 
ance included, is required, in order to better constrain the super- 
adiabatic efficiency of theoretical models. Moreover, we em- 
phasize that EBs are extremely useful tools for testing the valid- 
ity of stellar models, and consequently, precise data about such 
systems would be required. 



5 The case of VI 174 Ori is peculiar, and if the problem resides in 
the effective temperature determinations, then lithium abundances could 
also be affected by uncertainties greater than the quoted ones. 
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Table 2. Main parameters adopted for the selected set of EBs stars. 



System 


Mass [Mo] 


lo g r cff [K] 


log L/L Q 




[Fe/H] 


AS AS J05282 1+0338 5 (a) 


1 387 + 017 


3 708 + 009 


314 + 034 


3 10 + 20" 


-0.20 ± 0.20 e 


ASAS J05282 1+0338.5 (b) 


1.331 ±0.017 


3.663 ± 0.009 


0.107 ±0.034 


3.35 ±0.20" 


-0.10 ±0.20" 


EK Cep (a) 


2.020 ± 0.010 


3.954 ± 0.010 


1.170 ± 0.040 




+0.07 ± 0.05 f 


EK Cep (b) 


1.124 ±0.012 


3.755 ±0.015 


0.190 ±0.070 


3.11 ±0.30* 


+0.07 ± 0.05 f 


RXJ 0529.4+0041 A (a) 


1.270 ±0.010 


3.716 ±0.013 


0.140 ±0.080 


3.20 ± 0.30 c 


-0.01 ±0.04 * 


RXJ 0529.4+0041 A (b) 


0.930 ±0.010 


3.625 ± 0.015 


-0.280 ±0.150 


2.40 ± 0.50 c 


-0.01 ±0.04 * 


VI 174 Ori(a) 


1.009 ±0.015 


3.650 ±0.011 


-0.193 ±0.048 


3.08 ±0.20 d 


-0.01 ±0.04 s 


VI 174 Ori(b) 


0.731 ±0.008 


3.558 ±0.011 


-0.761 ± 0.058 


2.20 ±0.20'' 


-0.01 ±0.04 * 



Notes. For ASAS J052821.0338.5 (a) and (b) components we used the averaged value of [Fe/H] = -0.15 when computing the models. 

R eferences. (a) IStempels et all d2008h: {b) IMartin & Rebold d!993l) : {c) ICovino et al] d2000h : w iStassun et all d2004h : <e) IStempels et all d2008h ; 
^ iMartin & Rebold dl993h : M lD'Orazi et al] d2009fT ~ 



6. Conclusions 

We have discussed in detail the uncertainties on theoretical 
models by evaluating the effect of the errors affecting the ini- 
tial chemical composition and the up-to-date physical inputs 
(i.e. opacity, reaction rates, EOS), for several ages, masses, and 
chemical compositions. From these computations, we obtained 
a quantitative estimation of the error bar associated to the Li 
surface abundance predictions. The comparison between theory 
and observations was conducted on five young clusters, namely 
Ic 2602, a Per, Blanco 1, Pleiades, and Ngc 2516, and on four 
pre-MS EBs. 

Our results confirm the disagreement between present stand- 
ard models and 7 Li surface abundance in young stars. Motivated 
by the high sensitivity of 7 Li surface depletion to the convection 
efficiency and by the possibility of a dependency of the mixing 
length parameter on the evolutionary phase, gravity and T^, we 
therefore explored the effect of a different mixing length para- 
meter during the pre-MS and the MS evolution. We found that, in 
this case, very good agreement can be achieved for all clusters of 
the sample by adopting ffpMS = 1 -0, which is considerably lower 
than the MS one (i.e. ckms = 1-68 - 1.9), obtained by compar- 
ing theoretical isochrones with the observed colour-magnitude 
diagrams. 

We also checked the validity of such low-convection ef- 
ficiency models against four pre-MS detached double-lined 
eclipsing binaries, namely ASAS J05282 1+0338.5, EK Cep, 
RXJ 0529.4+0041 A, and VI 174 Ori. We found that, in the 
HR diagram, pre-MS tracks with a = 1.0 seem to agree well 
with the data for at least two of the three systems that have 
a star close to the Hayashi track (ASAS J05282 1+0338.5 and 
RXJ 0529.4+0041 A). However, the models are compatible with 
7 Li data both adopting the low- and high-convection efficiency, 
as a consequence of the large observational uncertainty still 
present on 7 Li abundance determinations. 

Our results point out the necessity of low-convection effi- 
ciency during the pre-MS phase in standard models. However, 
further analysis is required to clarify whether the mixing length 
parameter actually changes from low to high values evolving 
from the pre-MS to the MS phase, or if some physical mech- 
anism that acts in a way to partially inhibit con vection in the 
envel ope is lacking (i.e. magnetic field, see e.g., Ventur aet al j 
fl998h . 
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